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Abstract 
l> 

^* ■ We discuss the cosmological evolution of matter sources with small 

anisotropic pressures. This includes electric and magnetic fields, col- 
i | lisionless relativistic particles, gravitons, antisymmetric axion fields 

hip/ in low-energy string cosmologies, spatial curvature anisotropies, and 

stresses arising from simple topological defects. In many interesting 
cases the evolution displays a special critical behaviour created by 
Kn ' the non- linear evolution of the pressure and expansion anisotropies. 

C^ ■ The COBE microwave sky maps are used to place strong limits of 

order O a o < 5 x 10 -6 A (1 + z rec ) on the possible contribution 
of these matter sources to the total density of the universe, where 
1 < A < 3 characterises the anisotropic stress. In the case of cos- 
mological magnetic fields this translates into a present-day bound of 
Bq < 3.4 x 1CT 9 (f^o^Io) 5 Gauss. We explain why the limits obtained 
from primordial nucleosynthesis are generally weaker than those im- 
posed by the microwave background isotropy. The effect of inflation 
on all these stresses is also calculated. 



1 Introduction 

One of the key problems of cosmology is deciding which matter fields are 
present in the Universe. The most notable uncertainties govern the presence 
or absence of influential matter fields which might contribute to the total 
density of the Universe today, and thereby offer a resolution of the 'dark 
matter' [1] problem or provide an effective cosmological constant or vacuum 
stress [2] which alters the significance of the evidence governing the process 
of large-scale structure formation. Other matter fields, like magnetic fields 
or topological defects, can have a profound influence upon the evolution and 
properties of galaxies. In the early universe the identity of matter fields is 
even more uncertain. String theories are populated by large numbers of high- 
mass fields, but they will not survive to influence the late-time evolution of 
the Universe significantly if inflation takes place at high energies. Inflation, 
of course, involves it own material uncertainties: it requires the existence of 
a weakly-coupled scalar field which evolves slowly during the early universe 
[3]. In all these cases the principal problems are ones of fundamental physics 
- do the postulated matter fields exists or not?, what are the strengths of 
their couplings to themselves and to other fields? Relic densities residing in 
isotropic or scalar stresses are the simplest to determine. They can be cal- 
culated from a semi-analytic solution of the Boltzmann equation if their in- 
teraction strengths, lifetimes and masses are known. Limits on their masses, 
lifetimes, and interaction cross-sections are obtained by requiring that they 
do not contribute significantly more than the closure density to the Universe 
today, or produce too many decay photons in a particular portion of the 
electromagnetic spectrum [1], [4]. 

The most astrophysically interesting anisotropic stress is that of a cos- 
mological magnetic field. The origin of large-scale magnetic fields, whether 
observed in galaxies or galaxy clusters, is still a mystery. Intracluster fields 
are largely dominated by ejecta from galaxies. The invocation of protogalac- 
tic dynamos to explain the magnitude of the galactic field involves many 
uncertain assumptions but still requires a small primordial (pregalactic) seed 
field. Hence the possibility of a primordial field merits serious consideration. 
Other attempts to find an origin for the field in the early universe have ap- 
pealed to battery effects, the electroweak phase transition, or to fundamental 
changes in the nature of the electromagnetic interaction. All introduce fur- 
ther hypotheses about the early universe or the structure of the electroweak 
interaction. All aim to generate fields by causal processes when the Universe 



is of finite age. Therefore, any magnetic field created by these means will 
exist only on very small scales with an energy density that is a negligible 
fraction of the background equilibrium radiation energy density. 

Nevertheless, while such fields might still provide the seeds for non-linear 
dynamos in the post-recombination era, any large-scale magnetic field with 
a strength of order B ~ 10~ 8 Gauss, comparable to that inferred from the 
lowest measured intergalactic fields and close to the observational upper lim- 
its via Faraday rotation measurements, may well be of cosmological origin. 
A similar pregalactic (or protogalactic) field strength is inferred from the de- 
tection of fields of order 10~ 6 Gauss in high redshift galaxies and in damped 
Lyman alpha clouds, where the observed fields are likely to have been adia- 
batically amplified during protogalactic collapse. In the absence of a plausible 
dynamo for generating large-scale pregalactic fields, it is of interest to recon- 
sider the limits on a large-scale primordial field in view of new observational 
constraints that we outline below. 

Primordial magnetic fields can leave observable traces of their influence 
on the expansion dynamics of the Universe because they create anisotropic 
pressures and these pressures require an anisotropic gravitational field to 
support them. The influence of a magnetic field on reaction rates at nucle- 
osynthesis only limits the equivalent current epoch field to be less than about 
3 x 10~ 7 Gauss. We shall see that the cosmic microwave background isotropy 
provides a stronger limit on the strength of a homogeneous component of a 
primordial magnetic field. 

We will consider the general behaviour of anisotropic stresses, of which 
electric and magnetic fields are the simplest examples. We shall not as- 
sume that inflation has taken place, although we will also calculate the 
effects of inflation on their present abundances. The isotropy of the mi- 
crowave background requires any such anisotropies to be extremely small 
but we shall find that their evolution is subtle because of the coupling be- 
tween pressure anisotropies and the expansion dynamics during the radiation 
era. Anisotropic stresses have many possible sources - - besides cosmologi- 
cal magnetic and electric fields [5], examples are provided by populations of 
collisionless particles like gravitons [6], photons [7], or relativistic neutrinos 
[8], [9]; long-wavelength gravitational waves [10], [11], Yang-Mills fields [12], 
axion fields in low-energy string theory [13], [14]; and topological defects like 
cosmic strings and domain walls [15]-[17]. Our experience with high-energy 
physics theories also alerts us to the possibility that a (more) final theory 
of high-energy physics will contain many other matter fields, some of which 



may well exert anisotropic stresses in the universe. 

Cosmological observations of the isotropy of the microwave background 
tell us that the Universe is expanding almost isotropically. Therefore we can 
expect to employ the high isotropy of the microwave background to constrain 
the possible density of any relic anisotropic matter fields in the Universe. We 
shall see that a very general analysis is possible, largely independent of the 
specific identity of the physical field in question, which leads to strong limits 
on the existence of homogeneous matter fields with anisotropic pressures, 
irrespective of initial conditions in many cases. This analysis exhibits a 
number of interesting features of the general behaviour of small deviations 
from isotropy in expanding universes in the presence and absence of a period 
of inflation. 

In section 2 we set up the cosmological evolution equations for a universe 
containing a perfect fluid and a general form of anisotropic stress. Under 
the assumption that the anisotropy is small these equations can be solved. 
They reveal two distinctive forms of evolution in which the abundance of 
the anisotropic fluid is determined linearly or non-linearly, respectively. The 
evolution can be parametrised in terms of a single parameter which charac- 
terises the pressure anisotropy of the anisotropic fluid. In section 3 we give a 
number of specific examples of anisotropic fluids and explain how this analy- 
sis also allows us to understand the evolution of the most general anisotropic 
universes with 3-curvature anisotropy with respect to the simplest examples 
which possess isotropic 3-curvature. In section 4 we derive limits on the 
abundances of anisotropic stresses in the Universe today by using the COBE 
four-year anisotropy data [18] and describe the sampling procedure used to 
extract maximum information from the COBE maps. We also consider the 
case of an anisotropic stress which becomes non-relativistic after the end 
of the radiation era (analogous to a light neutrino species becoming effec- 
tively massive). In section 5 we consider the effect of a period of de Sitter 
or power-law inflation on the abundances of anisotropic stresses surviving 
the early universe, and in section 6 we compare the constraints we have ob- 
tained on anisotropic stresses from COBE with those that can be derived by 
considering their effects on the primordial nucleosynthesis of helium-4. The 
effect of the pressure anisotropy is to slow the decay of shear anisotropics to 
such an extent (logarithmic decay in time during the radiation era) that in 
the general case the COBE limits are far stronger than those provided by 
nucleosynthesis. 



2 The Cosmological Evolution of Skew Stresses 

Consider an anisotropic universe with metric [19] 

ds 2 = dt 2 - a 2 (t)dx 2 - b 2 {t)dy 2 - c 2 {t)dz 2 . (1) 

If we define the expansion rates by 

a = -; (3 = -; = -, 
a b c 

where overdot denotes d/dt, then the Einstein equations are 
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where T^ is the energy- momentum tensor and henceforth we set 8ttG = 1. 
We define the mean Hubble expansion rate, H, by 

H=^±, (6) 

and the two relative shear anisotropy parameters by 

R = — — and S = — — . (7) 

When R — S — the universe will be the isotropic flat Friedmann universe. 
We are interested in studying the behaviour of a non-interacting combi- 
nation of a perfect fluid and a fluid with an anisotropic pressure distribution. 
Thus the energy-momentum tensor for the universe is a sum of two parts: 

K = t a b + s a b , (8) 

where t% is the energy-momentum tensor of the perfect fluid, so 



t a b = diag(p^ -p m , -p m , -p„) (9) 

with equation of state 

Pm = ( 7 - ik, (10) 

and s b is the energy-momentum tensor of a fluid with density p and principal 
pressures pi,P2, and P3, so 

sl = diag(p,-p 1 ,-p2,-p3)- (H) 

We shall assume that the principal pressures of the anisotropic fluid are each 
proportional to its density, p, so that we have 

(Pl,P2,P3) = &P,VP,fJ>p), ( 12 ) 

where A, u, and p are constants. It will be useful to define the sum of the 
pressure parameters by 

A = A + /i + z/. (13) 

The ansatz of fli~2f) is clearly not the most general possible. We could 
seek to extend it by allowing A, u, and p to be time-dependent quantities. 
However, this opens up a huge range of complex, non-integrable behaviours. 
We know this because it allows the introduction of a Yang-Mills field which 
produces chaotic evolution (see Barrow and Levin [45]) equivalent to the 
motion of a point inside a potential bounded by four walls formed by the 
rectangular hyperbolae x 2 y 2 =constant. 

In the special case of an axisymmetric metric, (]l|), two of these three 
constants would be equal. The ratio of the densities of the two fluids is 
defined by 

Q = ~- (14) 

P* 

The two stresses t% and s£ are assumed to be separately conserved, so 

\7 b t a b = = V b s a b , (15) 



the a = components of which, using (li0|) -(|12|), give two conservation equa- 
tions, 



p, + 3#7P* = 0, (16) 

p + p[a(\ + l)+/3(ji+l)+9(v + l)] = 0. (17) 

Hence, the perfect fluid density falls in proportion to the comoving volume, 
as 

P * K jjj&> (18) 

and the anisotropic fluid density falls as 

1 

One can see from these two equations that Q oc a 1 ~ 1 b 1 ~ v ~ l c l ~ tl ~ l admits 
a variety of possible time-evolutions according to the values of the isotropic 
and anisotropic pressures. The evolution equations become, 

a + a 2 + a{(3 + 6) = -\p + Pl - p 2 - p 3 + (2 - j)p*} (20) 

i3 + P 2 + [3(a + 6) = i[p + P3-pi-P2 + (2- 7 )pJ (21) 

6 + 6 2 + 6(a + P) = I[p + p 2 -p3-p 1 + (2- 7 )pJ (22) 



Substituting H,R and S for a, (3, and 9, from (|7|) and using (pO|)-( p^) and 
(d), we have two propagation equations for the anisotropics, 

HR + RH + 3RH 2 = Pl - p 3 = (A - p)Qp„ (23) 

HS + SH + 3SH 2 = Pl - p 2 = (A - u)Qp„ (24) 

and the conservation equations flT6])-(|TTD combine to give 

Q = ^"[9(7 - 1) - 3(A + v + p) + R(2p -v-\) + S(2v -X-fi)] (25) 

The three equations (|23| ) -(|2~5| ) completely determine the evolution when the 
anisotropy is small. We are interested in the solutions of these equations in 



the case where the anisotropy is realistically small. Thus we assume that H 
and p^ take the values they would have in an isotropic Friedmann universe 
containing isotropic density p* (ie. S = R = Q = 0), so 



This is consistent with the ({]) — equation, (g). The three equations (|2^)-(P5D 
then completely determine the evolution of R, S, and Q. They reduce to 

fl + £ (2 _ 7) = ^zA, (27) 

S + 4(2-7) = ^>, (28) 

Q = 7^r[9(7-l)-3A + i?(2//-z/-A)+ (29) 
971 

+S{2u-X-fi)}. 



Axisymmetric solutions exist when S = and A = v. No essential simplifi- 
cation arises by imposing axial symmetry and so we shall treat the general 
case with S 7^ 0. 



The system of equations fl2~7|)-(|3~0"D gives rise to a characteristic pattern of 
cosmo logical evolution when stresses with anisotropic pressures are present. 
Equation (|30| ) reveals that there is a critical condition which, if satisfied, 
makes the problem a second-order stability problem; that is, if we linearised 
the equations about the isotropic solution (R — S — 0) with zero anisotropic 
stress density (Q = 0) we would find a zero eigenvalue associated with the 
evolution of Q(t). This critical condition is 

Criticality condition : 3(7 — 1) = A. (30) 

When 37 < 3 + A the shear distortion variables R and S relax towards their 
attractor where R — S — Oast — > 00, and so we have 

R = 2 « A -"> + 6 lt ^ -. 2Q } X ~»\ (31) 



2-7 2-7 

g(A-zz) 2=* 2Q(A- 

h dot t — > 

2-7 2-7 



2Q(A-z.) ^2 2Q(A - z/) 

5 = — - — — — + <5 2 t ^ -> — - — -— . (32) 



7-2 

with S\ and 62 constants. The Sit t terms are the contributions from the 
isotropic part of the 3-curvature. They fall off more rapidly than the part of 
the shear that is driven by the anisotropic pressure (assuming A 7^ p 7^ v). 
As t grows, the 5 terms become negligible if 37 > 2A (and recall that the 
isotropy is stable so long as 37 < 3 + A). Thus, R and S become proportional 
to Q(t), which is determined by 



= ^{9(7-l)-3A + ^-y[(A-/i)(2 / i-z,-A) + (A-z/)(2z,-A-/i)]}. 

(33) 
When the evolution is not critical (3(7 — 1) 7^ A) the evolution of Q(t) 
is determined by the terms linear in Q on the right-hand side of (|33"D; eqns. 
(|3l|)- (|32|) still hold, but now we have 



Q(t) = Q t k ; k^O, (34) 

k = A[3( 7 -l)-A]. (35) 

37 

Our assumption, eq.(p6p, that the evolution of p^ and H follow their values 
in an isotropic Friedmann universe will only be consistent at large times if 
k < 0, that is if 

-3< 3(7-1) < A< 3. (36) 

The right-hand side of this inequality derives from the causality conditions for 
signal propagation in the i = 1, 2, 3 directions (pi < p so A < l,/i< 1, v < 1); 
the left-hand side arises from pi > — p, which ensures the stability of the 
vacuum. When k > the anisotropic stress redshifts away more slowly 
that the isotropic perfect fluid on the average (over directions) and comes 
to dominate the expansion dynamics, making them completely anisotropic. 
We do not live in such a universe. By contrast, when k < 0, the isotropic 
stresses redshift away the slowest and increasingly dominate the dynamics, 
so the gravitational effect of the anisotropic stresses steadily diminishes. In 
the critical case the average stress energy of the anisotropic stress is counter- 
balanced by the isotropising effect of the perfect fluid and the shear evolution 
is determined by the second-order effects of the pressure anisotropy. This 
effect can be seen in the study of free neutrinos by Doroshkevich, Zeldovich 



and Novikov [8] and in the study of axisymmetric magnetic fields by Zeldovich 
[20]. 

When the evolution is critical, ( |30| ) holds, k = 0, and the evolution of 
Q{t) is decided at second-order in Q by 

AO 2 
Q = 9^ 2 )t « A -li)(2li-v-X) + {X- i/)(2i/ - A - n)}. (37) 

Hence, 

A = 9 7 (2 - 7 ) « A - Mfr - * - A) + (A - i/)(2i/ - A - //)}, (39) 

where Qo an d ^i are constants. For physically realistic stresses we have A > 0, 
and so as t — > oo the ratio of the energy densities approaches the attractor 

Q "> 77T7V (40) 

while the associated shear distortions approach the values given by eqns. 



(|31"1)-(|32|). Since the values of R and S at the epoch of last scattering of 
the microwave background radiation determine the observed temperature 
anisotropy we will be able to constrain the allowed value of Q at the present 
time by placing bounds on R and S at last scattering and then evolving the 
bounds forward to the present day. 

We can also calculate the asymptotic forms for the expansion scale factors 
of (|IJ). In the critical cases, as t — > oo, they evolve towards 



a{t) oc £^{ln£} m , (41) 

b(t) oc t^{\nt} w } (42) 

c(t) oc £^{ln£} n . (43) 

where the constants m, n, w are defined by 



m = 4(2A " " - V) (44) 

9 7 A(2- 7 ) ' l44j 
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= 4(2p- y -A) 

9 7 A(2- 7 ) ' v ; 



4 2z/-/i-A 
w = -5: — —T. ' 46 

9jA(2 - 7 v ; 



Some specific cases will be considered below. Note that the asymptotic form 
for the scale factors, fl4T|) -(f43D is consistent with the principal approximation 
imposed by eq. (EH). 

In the non-critical cases (with k < 0), as t — *■ oo, the scale factors evolve 
to first order as 

a(t) ex t^{l-V a t k }, (47) 

bit) oc t^{l-V b t k } } (48) 

c(t) oc t^{l-V c t k }, (49) 

where the constants V a ,Vb, and V^ are defined by, 

{V a , V b , V c } = 9fc ~^ x {/i + v - 2A, A + n - 2v, X + v - 2/i}. (50) 

These asymptotes reveal the different behaviour in the critical cases which 
produces the logarithmic decay of the shear. In the non-critical cases the 
anisotropic perturbations to the scale factors decay as power laws in time 
since k < 0. 



3 Some Particular Skew Stresses 

There are many examples of anisotropic stresses to which the analysis of the 
last section might be applied. We shall consider some of the most interesting. 
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In each case we can determine the characteristics of the 'critical' state in 
which the evolution of the anisotropy is determined by the non-linear coupling 
with the anisotropic stress. In section 4 we will go on to calculate the observed 
microwave background temperature anisotropy. The most familiar example 
of an anisotropic stress is that of an electromagnetic field. 

3.1 Magnetic or Electric fields 

In this case the anisotropic energy-momentum tensor s b a of eqn. (|ll|) has a 
simple form [20], [22], [24]. For example, for a pure magnetic field of strength 
B directed along the z-axis, we have Tq = T| = — T-j 1 = — T 2 2 = B 2 /87r and 
so this corresponds to the choice 

A = v = -ix = 1 (51) 

Hence, 

A(magnetic) = A(electric) = 1 (52) 

Therefore, from eqn.(|30D, we see that the criticality condition is obeyed 
when the background universe is radiation dominated. Thus, in the standard 
model of the early universe, the evolution of magnetic (or electric) fields 
will exhibit the non-linear logarithmic decay found in equations fl37D -fl40|). 



This was first pointed out by Zeldovich [20], and can be identified in the 
calculations of Shikin [21] and Collins [22]. Barrow, Ferreira, and Silk [23] 
used the COBE data set [18] to place constraints on the allowed strength of 
any cosmological magnetic field. 

As a specific example of a critical case, a radiation dominated universe 
(7 = 4/3) containing a magnetic field aligned along the z-axis, A = v — 
—\x = 1, gives A = 4 in (EM) and hence 



a(t) oc b(t) oc£5{l n £}4 (53) 

c{t) oc t?{\nt}-% (54) 

Notice that the volume expansion goes as in the isotropic radiation universe, 
abc oc £2 in accord with the principal approximation stipulated by (|26|). 

As an example of a non-critical case consider a pure magnetic field aligned 
along the z-axis of a dust universe (7 = 1). We have k = —2/3, and so at 
late times 
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2 40n 

a(t) oc b(t) oc t*{l - -^ + ...}, 

C(t) OC ta{l + ^ + ...}. 

This behaviour can be seen in the exact magnetic dust solutions of Thorne 
[24] and Doroshkevich [25]. Other studies of the evolution of cosmological 
magnetic fields in anisotropic universes can be found in refs. [28]. 

The axion field in low-energy string theory [13], [14] also creates stresses 
of this characteristic form and a source-free magnetic field has been used to 
study the possibility of dimensional reduction in cosmologies with additional 
spatial dimensions near the Planck time by Linde and Zelnikov [26] and 
Yearsley and Barrow [27]. 

3.2 General Trace-free stresses 

The magnetic field case is just one of a whole class of skew fields which exhibit 
non-linear evolution during the radiation era. Any anisotropic stress which 
has a trace-free energy-momentum tensor, s b a , will have 

A = 1 (55) 

and will exhibit critical evolution in the presence of isotropic black body radi- 
ation with 7 = 4/3. This case includes free-streaming gravitons produced at 
t P £ ~ 10~ 43 s which are collisionless at t > 10t P £ because of the weakness of the 
gravitational interaction mediating graviton scatterings [6]. It is also likely 
to include all asymptotically-free particles at energies exceeding ~ 10 15 GeV 
when interparticle scatterings, decays and inverse decays have interaction 
rates slower than the expansion rate of the universe, H. In all cases defined 
by A = 1 the shear to Hubble expansion rate decays only logarithmically 
during the radiation era, 

R oc S oc — t — r-. (56) 

In the dust era the critical condition will not continue to be met for trace- 
free fields. Although the evolution of Q, R, and S is now determined at linear 
order in (|3l|), (|32"D, and (|40|), there is still a significant slowing of the decay of 
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isotropy by the anisotropic stresses compared to the case where anisotropic 
stresses are absent. We see that when 7 = 1 we have 

— oc R oc S oc -—77 (57) 

compared to a decay of R oc S oc t^ 1 when anisotropic stresses are ab- 
sent (Q = 0) or ignored. The study by Maartens et al [29] of the evolution 
of anisotropics in a dust-dominated universe containing possible anisotropic 
stresses at second order (described by a tensor ir b a in reference [28] which is 
equivalent to the s b a used here) displays this same slowing of the shear decay 
(the discussion of ref. [41] omits this consideration). 

In general, we can see that if the criticality condition is satisfied in the 
radiation era, when 7 = 4/3, then it cannot be satisfied during the dust era 
that follows, when 7 = 1. However, an interesting situation can arise if the 
source of the trace-free stress is a population of particles which are relativistic 
above some energy (so A = 1 there) and non-relativistic when the universe 
cools below this energy scale (so A = there). Thus there can be a change 
in the value of A with time. We shall examine this case in section 4. 

Some matter fields with anisotropic pressures, like Yang-Mills fields [12], 
correspond to an anisotropic fluid with time- dependent A, /1, and v. But 
if the evolution is slow enough, it is well-approximated by the model with 
constant values of A, /1, and v used here. 

3.3 Long- wavelength gravitational waves 

In the last section we considered the evolution of anisotropic stresses in the 
simplest flat (f2 = 1) anisotropic cosmological model of Bianchi type I with 
the metric ([l]). The most general anisotropic universes of this curvature 
are of Bianchi type VIIq and they have an Einstein tensor that can be 
decomposed into a sum of two pieces: one corresponding to the Einstein 
tensor for the simple type I geometry, the other to a piece that describes 
additional long- wavelength gravitational waves [10], [11], [30]. These waves 
create anisotropics in the 3-curvature of the universe in addition to the sim- 
ple expansion-rate anisotropies present in the Bianchi I universe (which has 
isotropic 3-curvature). However, the contribution by the long- wavelength 
gravitational waves can be moved to the other side of the Einstein equa- 
tions and reinterpreted as an additional 'effective' energy-momentum tensor 
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describing a 'fluid' of gravitational waves. It has vanishing trace. This de- 
composition means that any general flat (or open) Bianchi type universe of 
type VIIq (or Vllh) containing an isotropic perfect fluid, ( |T0| ) behaves like 
a Bianchi type I (or type V) universe containing that fluid plus an addi- 
tional traceless anisotropic fluid. The parameters A, /i, and v are approxi- 
mately constant when the anisotropics are small. Hence the general evolution 
of anisotropic universes with anisotropic curvature containing isotropic ra- 
diation will exhibit the same characteristic logarithmic decay of the shear 
anisotropy given by eqns. fl3l"D , fl32|) , and ( flOf ) during the radiation era. 



This behaviour appears in Collins and Hawking [31], Doroshkevich et al [32], 
and other authors, [33], [34]. During the dust era any anisotropic curvature 
modes will behave like trace-free stresses and, although the evolution will no 
longer be critical, the shear anisotropy will fall more slowly (R oc S oc t~ 2 ' 3 ) 
than in simple isotropic universes with isotropic curvature and no anisotropic 
stresses (where R oc S oc t 1 ^ 27 ). 

3.4 Strings and Walls 

Topological defects provide another class of anisotropic energy-momentum 
tensors which fall into the category of stress modelled by (|i~2"|) for part of their 
evolution. The energy-momentum tensors of string sources were also consid- 
ered more generally by Stachel [16], Marder and Israel [17]. The specific de- 
scription of line stresses created by topological defects is reviewed in ref. [15]. 
An infinite string with mass per unit length [i extending in the x-direction 
contributes an anisotropic stress-tensor s b a = fi5(z)5(y)diag(l, 1, 0, 0); that 
is 

A = -1, n = v = -> A(string) = -1. (58) 

Their evolution could therefore only be critical in a universe containing 
a perfect fluid with equation of state p* = — p*/3. It is worth noting that 
this corresponds to the evolution of the curvature term in the Friedmann 
equation when the universe is open. Thus we would expect the string stress 
to evolve critically during the late curvature-dominated stage (1 + z < Qq 1 ) 
of an open universe, during a curvature-dominated pre-inflationary phase, or 
during any period when quantum matter fields with 7 = 1/3 dominate the 
expansion. 

For slow-moving infinite planar domain walls of constant surface density 
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in the x — y plane , the stress corresponds s b a oc r)(z)diag(l, 1, 1, 0), where 
rj(z) is a local bell-shaped curve centered on z = [15], and corresponds to 
the choice 

A = v = -1, /U = -»■ A(waZZ) = -2. (59) 

The wall stress would only be critical if the equation of state of the back- 
ground universe is p* = — 2p^/3. 

These descriptions do not include the complicated effects of the non- 
linear evolution of a population of open strings and loops which must include 
intersections, kinks, gravitational collapse and gravitational radiation. 

3.5 General classification and energy conditions 

When evaluating the form of the residual anisotropy and energy density 
created by anisotropic cosmological stresses it is most convenient to classify 
stresses by the value of A. We can circumscribe the likely range of realistic 
A values by considering some general restrictions on the energy-momentum 
tensor. The dominant energy conditions [35] require us to impose the physical 
limits 

|A| < 1, \ft\ < 1, \u\ < 1. (60) 

Hence we have the bounds 

- 3 < A < 3. (61) 

If the strong-energy condition [35] were imposed on the energy-momentum 
tensor s b a then we would have a stronger restriction 

A > -1. (62) 

This condition is violated by string and wall stresses and necessarily by any 
isotropic stress which drives inflation (0 < 7 < 2/3) since A(wall) = —2 
and A(string) = — 1. We see from fl34|)-(|35D that there can only be approach 
to isotropy at late times in the radiation era if A > 1. When < A < 1 
the expansion approaches isotropy during the dust era but not during the 
radiation era. When A < isotropy is unstable during both the dust and 
radiation eras. 
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4 Microwave background limits 

Let us consider the simplest realistic case in which the Universe is radiation 
dominated until l + z eq = 2.4 x 10 4 f2o/i(j an< ^ then dust dominated thereafter; 
here; Qq < 1 is the cosmological density parameter and ho is the present 
value of the Hubble constant in units of lOOKms' 1 Mpc~ l . We shall assume 
that the microwave background was last scattered at a redshift z rec , where, 
in the absence of reheating of the cosmic medium, 

1 + z rec = 1100. (63) 

If there is reionization of the universe then last scattering can be delayed 
until 1 + z rec = 39(fi w /i ) _1 for Qoz rec « 1, where f2 b0 is the present baryon 
density parameter [36]. 

The microwave background temperature anisotropy is determined by the 
values of shear distortions R and S at the redshift z rec . The evolution of the 
photon temperature in the x, y, and z directions is given by 



T x = T —^-=T exp{-adt}, (64) 

a{t) J 

T y = T -^ = T exp{-Jpdt}, (65) 

T z = T o ^- = T o exp{-j0dt}. (66) 

If we define the temperature anisotropy by 

8T _{T X -T y ) +{T X -T Z ) 



, T (67) 

then, for small anisotropies (so exp{— J adt} ~ 1 — J adt} etc), using the 
definitions of (0) in (|6"TD, we have 

5T r 

— = -H J(R + S)dt (68) 

Since recombination always occurs in the dust era, H = 2/3t,and the ob- 
served microwave background anisotropy will be 



ST 
T 



R + S 

A , 



xm<f>,Sl ). (69) 
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where / ~ 0(1) is a pattern factor taking into account non-gaussian statis- 
tical factors and the possible pattern structure created by complicated forms 
of anisotropy in the general case [13], [23], [37]. Typically, in the most gen- 
eral homogeneous flat universes the pattern combines a distorted quadrupole 
with a spiral geodesic motion. In ref. [23] a detailed discussion of the sam- 
pling statistics of the microwave background temperature distribution and 
the non-gaussian nature of the anisotropic pattern was given. This enables 
limits to be derived from the whole COBE sky map rather than, say, simply 
from the quadrupole as is generally done. 

The problem of constraining global anisotropy is substantially different 
from the traditional statistical task of estimating parameters in Gaussian 
models. In the latter case, the ensemble is entirely characterized by the 
power spectrum while in the former, a given set of parameters corresponds 
to a completely specified pattern in the sky, up to an arbitrary rotation. This 
problem was dealt with in some detail in Bunn et al [33] . A brief outline of 
our procedure is as follows. One can model the microwave background signal 
as the sum of two components: a statistically isotropic Gaussian random 
field AT/, which we assumed to have a scale invariant power spectrum on 
the scales we are interested in, and a global, anisotropic pattern, ATa, which 
is uniquely defined by the set of parameters x (which measures the spiral 3- 
curvature anisotropy), Qq, /150, (a/H) , and 6, (its orientation on the sky). 
Each pixel of a data set of measured microwave background anisotropics is 
given by d,i = (AT*/3)(rj) + iVj where /3 represents the DMR beam pattern, 
Yi is a unit vector pointing in the direction of pixel i, iVj is the noise in pixel 
% and * is the convolution operator. To an excellent approximation, one can 
assume that N is Gaussian "white" noise, i.e. (NiNj) = s^Sy. 

Our task is, given a pair (x, Qq), to find the orientation (6, 0) which al- 
lows the maximum observed value of (o~/H) . One can do this using standard 
frequentist statistical methods: we define a goodness-of-fit statistic that de- 
pends on the data, compute its value for the actual data, and then compute 
the probability that a random data set would have given a value as good as 
the actual data. In Bunn et al [33], rj was defined to be: 



i] = min^; where r] 1 = — ^— -g — -; (70) 



A 2 -A 2 
M^'' 1 ' '"■ A 2 

Aq is the noise-weighted mean-square value of the data and A 2 is the noise- 
weighted mean-square value of the residuals after we have subtracted off the 



anisotropic part. Note that removing the incorrect anisotropic portion will 
only increase the residuals so the difference between the two terms is an ob- 
vious choice for a goodness-of-fit. Dividing by Aq ensures a weak dependence 
on the amplitude of the isotropic component, while defining the statistic as 
the minimum of r] l allows us to deal with the uncertainty in (6, </>). 

This statistical method was applied to the 4-yr COBE DMR data-set. 
The two 53 GHz and the two 90 GHz maps were averaged together, each 
pixel weighted by the inverse square of the noise level, to reduce the noise 
level in the average map. All pixels within the Galactic cut were removed so 
as to reduce Galactic contamination, and a best-fit monopole and dipole were 
subtracted out. The map was degraded from pixelization 6 to pixelization 5 
(i.e. binning pixels in groups of four). Simulations were performed for a set 
of models from the (Qq, x) plane; for each choice of the three parameters (Qq, 
x and (a/H)o) approximately 200 to 500 random DMR sets were generated, 
so allowing us to determine an approximate fit to the probability distribution 
function of rj. 

This leads to bounds on / in flat and open universes of 

0.6</<2.2. (71) 

This bound can be improved by a factor of \/3 if one considers the results 
from Kogut et al [42] . In this case, a slightly different goodness-of-fit statistic 
is used: instead of working with the noise- weighted quantities, Aq and Af, 
the authors chose to weight the pixels with the covariance matrix of the total 
Gaussian components (i.e. the noise and isotropic cosmological components). 

The discussion above shows that the A > 1 case is the realistic one which 
allows evolution towards isotropy at late times. The observed anisotropy is 
therefore given in terms of the present value of the density ratio, Q(t ), by 

5T _ 2Q(t )f(2\-y-v)(l+z rec ) A 
T ~ A {fZ) 

If we take the COBE 4-year data set to impose a limit of ST/T < 10~ 5 on 
contributions by the anisotropic fluid, and use (|7T|), then the present density 
parameter of the anisotropic fluid, Q a0 , is limited by 



< 8.3 x 10- 6 A = / A \ / 1100 \ A 10 _ 5 .08-3.04A 

a0 - (2\- fi-u)(l + z rec ) A \2X- fi-u) \ l + z rec J 

(73) 
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Since 1 < A < 3 for realistic fluids we can examine the extremes of 
this limit which is strongest when there is no reionization of the Universe at 
z « 1100 because reionization allows a longer period of power-law decay 
of R, S, and Q, hence a smaller residual effect on the microwave background 
isotropy. In the two extreme cases we have limits of 

A = 1 : n a0 < 3.7 x KT 9 x (^-) (74) 

V 1 T" Z rec ' 

/ 1 1 00 \ 3 

A = 3 : n a0 < 9.5 x 10~ 15 x (75) 

V i + z rec / 

The A = 1 case with the choice of (^I|) corresponds to limits on the present 
cosmological energy density allowed in magnetic (or electric fields) which was 
studied in [20] and in [23] where the most general evolution of anisotropy was 
included. Note that these limits are far stronger than those that are generally 
obtained for isotropic forms of dark matter by imposition of astronomical lim- 
its on the maximum total density, the age of the Universe, or the deceleration 
parameter [1], [4]. 

The most conservative limit on the cosmological magnetic field arises 
when we assume that the whole anisotropy is contributed by the magnetic 
field stresses. This gives us a final bound on the magnitude of the magnetic 
field today of This gives a limit of 

B < 2.3 x 10- g f^ (O /4)5Gauss (76) 

Adams et al [43] have used the weak nucleosynthesis limit of Grasso and 
Rubinstein of B < 3 x 10~ 7 Gauss obtained for random fields to argue 
that a cosmological magnetic field could lead to observable distortions of the 
acoustic peaks in the microwave background. Our limit on B Q rules out any 
observable effect of a homogeneous magnetic field on the acoustic peaks. A 
large scale, inhomogeneous, magnetic field may, however, introduce observ- 
able distortions in the acoustic peaks. Our limit permits a field strength of 
10 -9 Gauss required to induce a measurable Faraday rotation in the polar- 
ization of the microwave background [44]. 

The limit ([76]) is not strong enough to exclude a cosmological origin for 
galactic magnetic fields from purely adiabatic compression of a background 
field down to galactic scales with no significant enhancement from a dynamo. 
If that were the source of galactic fields then one would expect significant 
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fields to exist in elliptical galaxies as well as spirals. If galactic fields were 
amplified significantly by dynamo action then one would not expect signifi- 
cant ordered fields to exist over large scales within ellipticals, although they 
could exist in local disordered forms following ejection from stars. 

4.1 Evolution of anisotropic dark matter with a char- 
acteristic energy scale 

We are familiar with the standard picture of the cosmological evolution of 
light (m << lMeV) weakly interacting particles [1], [4]. When T > m 
they behave like massless particles with their number density similar to that 
in photons, up to statistical weight factors. Their energy density redshifts 
away like (1 + z) 4 until the temperature falls to the value of their rest mass. 
Thereafter they behave like non-relativistic massive particles and their en- 
ergy density redshifts away more slowly, as (1 + z) 3 . We can consider an 
analogous scenario for anisotropic stresses. Suppose we have an anisotropic 
fluid which behaves relativistically until the temperature falls to some value 
T + < lO^K m leV and then behaves non-relativistically at lower tempera- 
tures. This situation corresponds to following the evolution with A = 1 for 
T > T + and then with A = for T < T + . The evolution is complicated 
by the fact that it will be critical during the radiation era but non-critical 
during the first stage of the dust era when the temperature exceeds T + when 
Q will decay in accord with (0), before becoming critical again during the 
remainder of the dust era until the present when Q will decay logarithmi- 
cally in accord with fl4"0|). In this case the final density of the anisotropic 
dark matter is constrained by the microwave background isotropy to have 



where 



n a0 < 1.2/- 1 x io- 4 (-^2L) x - -L (77) 



A = -{(// - X)(2/j -u-X) + (u- X)(2v - A - fj,)} (78) 

9 

For example, if the anisotropic matter has a characteristic energy scale m 
then, since m = T + = 2.4 x 10 _4 (l+2; + ) eV, we have (picking A = v = —pL = 1 
so A = 32/9), with (171]), that, 
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^ao < 5.7 x lO" 5 (j^-) x \ (79) 

\l + z rec J (8.3 + In M e y) 

and there is a very weak dependence on the mass scale m e y = (m/leV). For 
m > leV we have, roughly, that 

n a0 < 6.9 x 10- 6 (t^-) (80) 

Therefore these fields are always constrained by the microwave back- 
ground anisotropy to be a negligible contributor to the total density of the 
Universe. 



5 The Effects of Inflation 

So far we have ignored the consequences of any period of inflation in the very 
early universe [3]. The equations (|27|)-(|3TD hold for the case of generalised 
(power- law) inflation with < 7 < 2/3, [37]. However, for de Sitter inflation 
with 7 = 0, the isotropic density drives the inflation with p^ = 3Hq = 
constant, and they are changed to 



R + 3RH = 3# Q(A-/i), (81) 

S + 3SH = 3H Q Q(X-u), (82) 

Q = -Q^-[-9-3A + R(2v-v-\)+S(2v-\-v)].(83) 

This system can only be critical for the evolution of Q only if A = —3 
but this cannot occur for an anisotropic fluid subject to (|60|). As t — > 00 we 
have 



R -► q{\ - //) + 5x exp(-3H t) (84) 

S -»■ q(\-v) + 5 2 exp(-3H t) (85) 

Q -> g exp[-(3 + A)/Jot] (86) 

Thus, since A + 3 > the contribution made to the shear by the anisotropic 
pressure decays. Moreover, we see that the contribution by the isotropic 
curvature (5) terms to the shear dominates the contribution by the pressure 
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anisotropy at late times if A > 0. In both cases the shear anisotropy decays 
away exponentially fast. This is in accord with the expectations of the cos- 
mic no hair theorems [38] because the anisotropic stresses considered here 
obey the strong energy condition when A > —3. Thus if N e-folds of de 
Sitter inflation occur in the very early universe, then the values of the shear 
distortion parameters, R and S, are each depleted by a factor exp(— 3N) 
whilst the relative density in the anisotropic fluid, q, is reduced by a factor 
exp{-(3 + A)iV}. 

Similarly, in the case of power-law inflation (0 < 7 < 2/3), the isotropic 
curvature mode will dominate the late time evolution of the shear during the 
inflationary phase if 



2A , N 

7 < -3-. (87) 



If anisotropic stresses can be generated at the end of inflation then the 
analysis of the previous sections applies. 

6 Primordial Nucleosynthesis 

It is instructive to consider the question of whether primordial nucleosyn- 
thesis can provide stronger limits than microwave background isotropy on 
the densities of anisotropic stresses in the Universe. The key issue that de- 
cides this is whether the microwave background constraint on the expansion 
anisotropy created by an anisotropic energy density, which is of order 10 -5 
and imposed at a redshift z rec ~ 10 3 or lower (in the event of reionization), 
is stronger than a limit of order 1 — 10 _1 on changes to the mean expansion 
rate imposed at the epoch of neutron-proton freeze-out, Zf r ~ 10 10 . The is- 
sue is decided by knowing how fast the anisotropic stresses decay with time 
between Zf r and z rec . 

In the simplest anisotropic universes, which have isotropic 3-curvature, 



the shear anisotropy falls like the that of the 5 mode in eqns. (|3T|) -(|32l). 
Consider first the simple isotropic curvature case with no anisotropic stress, 
so Q = 0. The evolution of the shear to Hubble rate parameters, R and S, 
during the dust and radiation era of a universe with Q = 1 is 
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z > z eq : R oc S oc t 2 oc 1 + z (88) 

z < z eq : RocSocT 1 oc (l + z)i (89) 

Hence, if nucleosynthesis gives an upper limit of ~ 0.2 (roughly equiv- 
alent to adding one neutrino type to the standard three-neutrino model) 
on the values of \R\ and \S\ at the redshift of neutron-proton freeze-out, 
Zf r , this corresponds to an upper limit at the time of last scattering of 
0.2(1 + Zf r )~ l {\ + z eq ) 2 (1 + z rec )2, so we have 

3 

R rec ~ S rec < 4.6 x 10" 9 (^o)" x (-j^p) ' (90) 

Since the microwave background gives a limit of ST/T ~ (R + S% ec < 
10 _5 ,[18], we see that the nucleosynthesis limit is much stronger when the 
3-curvature is isotropic, as first pointed out by Barrow [39]; nor are conceiv- 
able improvements in microwave receiver sensitivity ever likely to close to gap 
of 10 4 that exists between the nucleosynthesis and microwave limits on the 
isotropic 3-curvature 5 modes. These effects of these simple anisotropy modes 
on nucleosynthesis were considered in the papers of Hawking and Tayler [40] 
and Thorne [24]. 

However, the situation changes when we consider the most general forms 
of anisotropy with anisotropic curvature or when matter is present with 
anisotropic pressures. If the evolution will be critical during the radiation 
era the anisotropy falls off so slowly during the period before z eq that the 
microwave background limits become stronger than the nucleosynthesis lim- 
its. This effect was considered in both the dust and radiation eras by Barrow 
[39] in the evolution of Bianchi type VII universes close to isotropy but the 
observational limits were 100 times weaker then. For example, in the in- 
teresting case where A = 1 the evolution of anisotropies follows the form, 



z > z eq : R oc S oc yj£ y (91) 

z < z eq : R oc S oc 1 + z. (92) 

Hence, the nucleosynthesis limit translates into an upper limit on R and S 
at the redshift of last scattering of only 
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Rrec~S rec <7.lxlO- A (n h 2 )- l x (^y^p)- (93) 

We see that this is never stronger than the microwave background limit of 
< 10~ 5 for any possible redshift of last scattering. Alternatively, we might 
restate this result as follows: it is possible for anisotropic fluids to create a 
measurable temperature anisotropy in the microwave background radiation 
without having any significant effect upon the primordial nucleosynthesis of 
helium-4. 



7 Discussion 

We have described the cosmological evolution of a very general class of 
anisotropic stresses in the presence of an isotropic perfect fluid. Such anisotropic 
stresses encompass a wide range of physically interesting cases, including 
those of cosmological electric and magnetic fields, a variety of topological 
defects, gravitons, and other populations of collisionless particles. When 
the universe is close to isotropy the mean expansion rate behaves as in 
the isotropic Friedmann universe to leading order, but the density of the 
anisotropic stresses is coupled to the expansion anisotropy in an interesting 
way. It was found that there are two possible forms for the evolution. In 
the critical case the density of anisotropic matter is determined by its non- 
linear coupling to the expansion anisotropy and density of anisotropic stress 
falls only logarithmically relative to the isotropic background density. By 
contrast, in the non- critical case the density falls as a power-law in time 
relative to the background density. In all cases the microwave background 
anisotropy can be used to place a limit on the density of matter that could 
be residing in the universe today in forms with anisotropic pressures. These 
limits arise because their evolution approaches an asymptotic attractor in 
which the temperature anisotropy produced in the microwave background 
by the anisotropic stresses is simply related to their density relative to the 
isotropic background density. The limits obtained on the possible cosmolog- 
ical density of matter of this form are far stronger than conventional limits 
on isotropic forms of dark matter because they make use of the microwave 
isotropy limits (5T/T < 10~ 5 ) rather than the far weaker limits from obser- 
vations of the Hubble flow and age of the Universe (f2 < 0(1)). We went on 
to explain why the limits that can be impose upon anisotropic stresses by the 
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microwave background isotropy measurements are in general stronger than 
those arising from the limits on the primordial nucleosynthesis of helium-4. 
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